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A new universal empirical function that depends on a single critical exponent (acceleration ex- 
ponent) is proposed to describe the scaling behavior in a dissipative kicked rotator. The scaling 
formalism is used to describe two regimes of dissipation: (i) strong dissipation and (ii) weak dissipa- 
tion. For case (i) the model exhibits a route to chaos known as period doubling and the Feigenbaum 
constant along the bifurcations is obtained. When weak dissipation is considered the average action 
as well as its standard deviation are described using scaling arguments with critical exponents. The 
universal empirical function describes remarkably well a phase transition from limited to unlimited 
growth of the average action. 
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I. INTRODUCTION 

In 1969 Boris Chirikov [l|, Q proposed what became 
one of the most important and extensively studied sys- 
tems in nonlinear dynamics and in the theory of Hamil- 
tonian systems and area-preserving maps [sl, Q , namely 
the Chirikov standard map. The model describes the 
motion of a kicked rotator. Applications of the model 
can be made in different fields of science including solid 
state physics |5|] , statistical mechanics [6] and accelerator 
physics It has also been studied in relation to prob- 
lems of quantum mechanics and quantum chaos [8|, Q, 
plasma physics [lo| and many others. 

The standard map is a dynamical system in which 
the nonlinearity given by a sine function is controlled 
by a parameter K. In the absence of dissipation, if K is 
small enough the structure of the phase space is mixed in 
the sense that Kolmogorov-Arnold-Moser (KAM) invari- 
ant tori and islands are observed coexisting with chaotic 
seas pj]-[l5| . As the parameter K increases and becomes 
larger than Kc ~ 0.971635 • • the last invariant span- 
ning curve disappears and the system presents a global 
chaotic component where a chaotic orbit spreads over the 
phase space. However, the introduction of dissipation in 
the model changes the mixed structure and the system 
exhibits attractors (l6l4l9j. We consider two different 
ranges of dissipation namely: (i) strong dissipation - the 
situation where the variable action loses more than 70% 
of its value upon a kick - where period doubling bifurca- 
tion cascade is observed and the Feigenbaum coefficient 
is numerically obtained and; (ii) weak dissipation where 
the average action exhibits scaling features [20|. Here 
we propose a new universal function that describes the 
scaling behavior of the average action by the knowledge 
of a single critical exponent, the so called acceleration 
exponent. 

The paper is organized as follows. In section |II] we 
describe all the necessary details to obtain the two- 



dimensional map that describes the dynamics of the sys- 
tem and also we present and discuss our numerical re- 
sults. Conclusions are drawn in Sec. IIIII 



II. THE MODEL AND RESULTS 

The Hamiltonian that describes the standard map has 
the following form (see e.g., [l, S HH ) : 



t2 oo . 

H=- + Kcos{0) J2 ^[f 



(1) 



where K is the amplitude of the delta-function kicks 
(pulses) with period T = 27t/u. The equation of motion 
is given by 



(2) 



Assuming that {In-, On) be the values of the variable im- 
mediately after the n^^ kick, (/n+i,^n+i) represent their 
values after the (n + 1)^^ kick, and introducing the dis- 
sipation parameter 7 the dissipative standard map 
is written as 



(1 



Ksm^On) 

mod{27r) 



(3) 



where 7 G [0, 1] is the dissipation parameter. If 7 = 
all the results for the Hamiltonian area-preserving stan- 
dard map are recovered. On the other hand, if 7 = 1 
the results for the one-dimensional sine-circle map are 
obtained [3]. In this paper we shall consider < 7 < 1. 
The system is area preserving only when 7 = since the 
determinant of the Jacobian matrix is det(J) = 1 — 7. 
Another consequence of the dissipation is that the map 
(j3j) is periodic only in 6>, while the action / ranges be- 
tween —00 < / < 00. 
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FIG. 1: Bifurcation cascade for (a) / and (b) both plotted 
against the parameter K. In (c) it is shown the Lyapunov 
exponent associated to (a) and (b). The damping coefficients 
used for the construction of the figures (a), (b) and (c) were 
7 0.80. 



FIG. 2: (Color online) (a) Phase space for the conservative 
standard map (7 = 0) with K — 1. (b) Basin of attraction for 
the attracting fixed points (sinks) of period 1 (red and black), 
2 (cyan), 3 (maroon) and 4 (green). The control parameters 
used to construct the basin of attraction were K = 1 and 
7= 10-^ 



In this paper we will consider two situations for 7, 
namely, strong dissipation and weak dissipation. The 
strong dissipation refers to the situation in which the ac- 
tion variable / loses more than 70% of its value upon a 
kick. We have considered the case 7 = 0.80. To explore 
some typical behavior we have used as initial conditions 
^0 = 3 and /q = 5.53 and investigated its attraction to 
periodic orbits, and looked at the bifurcations as K varies 
in the range where we have global chaos if 7 = 0. Fig. 
[TJa) shows the behavior of the asymptotic action plot- 
ted against the control parameter where a sequence 
of period doubling bifurcations is evident. A similar se- 
quence is also observed for the asymptotic variable 6>, as 
it is shown in Fig. [IJb). There are three small sequences 
of bifurcations, however we are not interested in their be- 
havior here, but the same procedure can also be applied 
to them. Observe that, of course, the bifurcations of the 
same period in (a) and (b) happen for the same value of 
the control parameter K. As one can see in Fig. [1] (c), 
as the parameter K approaches the point of bifurcation. 



the Lyapunov exponent approaches zero. As discussed in 
(23! the Lyapunov exponents are defined as 



A. 



lim — In |Aj 



J = 1,2 



(4) 



where Kj are the eigenvalues of M = fir^i ^ii^i^^i) 
Ji is the Jacobian matrix evaluated over the orbit (J^, Oi). 
If at least one of the Xj is positive then the orbit is clas- 
sified as chaotic. We can see in Fig. ([1]) (c) the behavior 
of the Lyapunov exponents corresponding to both Figs. 
([l])(a,b). It is also easy to see that when the bifurca- 
tions happen, the exponent A vanishes. Such a behavior 
occurs because the eigenvalues of the Jacobian matrix 
become complex numbers on unit circle. The Lyapunov 
exponents between 5.66 < K < 5.69 correspond to the 
small sequence of bifurcation observed in Fig. [1] (a,b) for 
the same range of the control parameter K. It is im- 
portant to emphasize that in some regions the Lyapunov 
exponent assumes a constant and negative value. 
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The period doubling cascade was discovered by May 
[23 in 1976 and by Grossmann and Thomae in 1977 
[25] and later it was discovered by Feigenbaum [26*, ^2f\ 
that there is a universal feature along the bifurcations. 
The period doubling bifurcations converge geometrically 
to the chaos border. The procedure used to obtain the 
Feigenbaum constant 6 is as follows: let Ki represent the 
control parameter value at which period- 1 gives birth to a 
period-2 orbit, K2 is the value where period-2 changes to 
period-4 and so on. In general the parameter Kn corre- 
sponds to the control parameter value at which a period- 
2"^ orbit is born. Thus, the Feigenbaum's 5 is written 
as 



. T — Kn-l 

= lim — 



(5) 



The theoretical value for the Feigenbaum constant S is 
S = 4.669201609 • • • . Considering the numerical data 
obtained through the Lyapunov exponents calculation, 
the Feigenbaum's 6 obtained for the dissipative standard 
map is (5 = 4.6691(1). Since the numerical results are 
very hard to be obtained for bifurcations of higher or- 
ders, we have considered in our simulations only the bi- 



furcations up to eleventh order, but still our result is in 
a good agreement with the Feigenbaum's universal 6. 

As we have shown for strong dissipation the system 
exhibits a sequence of period doubling bifurcations and 
for one of them we have found the Feigenbaum's S. Now, 
consider the parameter K such that the conservative dy- 
namics has the mixed phase space. When dissipation is 
taken into account the mixed structure is changed. Then, 
an elliptical fixed point (generally surrounded by KAM 
islands) turns into a sink. Regions of the chaotic sea 
might be replaced by chaotic attractors. In Fig. [2fa) the 
behavior of the phase space for the conservative dynam- 
ics (7 = 0) with = 1 is shown. As it is well known 
for such a value of the parameter K the last invariant 
torus is destroyed and the phase space has only KAM 
islands and a large chaotic sea. Figure Efb) shows the 
basin of attraction where the main fixed points are of 
period 1 (red and black), 2 (cyan), 3 (maroon) and 4 
(green). The dissipation parameter is 7 = 10~^. The 
procedure used to construct the basin of attraction was 
to divide both / and into grids of 1000 parts each, thus 
leading to a total of 10^ different initial conditions. Each 
initial condition was iterated up to n = 5 x 10^. Apart 
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FIG. 5: (Color online) (a) Behavior of cj(n,7) as a function 
of n for different values of the control parameter 7. (b) Their 
collapse onto a single universal plot. 



from the mentioned periodic attractors there are many 
others, which, however, are difficult to find and display, 
because their basins of attraction are very small. A more 
systematic overview of the dynamical phase diagram in 
the parameter-space (i^, 7), is given in reference 



When the control parameter K is sufficiently large, all 
the regular regions are destroyed and the phase space for 
the conservative case is fully chaotic. In this sense we 
will study some dynamical properties in the regime of 
K » Kc but taking into account weak dissipation. Our 
numerical results concern basically the behavior of the P 
which is also the average energy {E = P /2). Two steps 
were applied in order to obtain P. Firstly, we evaluate 
P over the orbit for a single initial condition 



j=0 



(6) 



where the index i corresponds to a sample of an ensemble 
of initial conditions. Hence P is written as 
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FIG. 6: (Color online) (a) Plot of Usat as a function of 7. (b) 
Behavior of the crossover number Ux vs. 7. 



We have iterated Eq. ([7]) up to n = 10^ for an ensemble 
of M = 2000 different initial conditions. The variable Iq 
is fixed at Iq = 10~^i^^, ^0 is uniformly distributed on 
^0 ^ [0, 27r] and the control parameter was fixed as = 
20. Figure [3] shows the behavior of the P as a function of 
n for five different values of the dissipation parameter 7, 
as labelled in the figure. It is easy to see in Fig. [3] (a) two 
different kinds of behavior. For short n, P(n^^) grows 
according to a power law and suddenly it bends towards 
a regime of saturation for large enough values of n. The 
crossover from growth to the saturation is marked by a 
crossover iteration number Ux . This crossover number Ux 
is actually quite well defined by the crossing of the two 
straight lines: The acceleration line and the saturation 
line in the log- log plot. Based on the behavior shown in 
Fig. [3] (a), the following scaling hypotheses are assumed: 

1. When n <C n^c, /^(n,7) grows according to 

P(n,-i)(xn^ , (8) 
where the exponent p is the acceleration exponent; 

2. When P{n^ 7) approaches a regime of sat- 
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uration, which is described by 

P satin, , (9) 

where the exponent a is the saturation exponent; 

3. The crossover iteration number that marks the 
change from growth to the saturation is written as 



10^ 



Tlx OCJ 



where z is the crossover exponent. 



(10) 



After considering these three initial suppositions, 

P{n,j) can be described in terms of a scahng function 
of the type 



where r is the scahng factor, p and q are scahng expo- 
nents that in principle must be related to a, P and z. 
Since r is a scaling factor we can choose r^n = 1, or 
r = n~^^P. Thus, Eq. (pT]) is rewritten as 



(12) 



where the function Pi(n~^^^^) = P{l,n~^^^j) is as- 
sumed to be constant for n <^ Ux- Comparing Eqs. ([8j) 
and ([12]) , allows one to conclude that {3 = — Choos- 
ing now r = 7~^/^ and Eq. (pTj) is given by 



/2(n, 7) = 7-^/^/2^(7-^/%) , 



(13) 



where p2{j~^^^n) = /2(7~^/^n, 1). It is also assumed 
as constant for n ^ Ux- Comparison of Eq. ([9]) and 
Eq. (p!3|) gives us a = —l/q. Given the two different 
expressions of the scaling factor r in the crossover region 
n ~ n^;, we have r = = 7" and one must conclude 
that the crossover exponent z is given by z = a/ j3. Ob- 
serve that the scaling exponents are determined if the 
critical exponents P and a were numerically obtained. 
The exponent j3 is obtained from a power law fitting for 
P{n,j) curves for the parameter 7 G [10-^,10-^] for 
short iteration number, n « Ux. Thus, the average of 
these values give us /3 = 1.004(3). Figure |4] shows the be- 
havior of (a) P satvs.^ and (b) UxVs.j. Applying power 
law fitting in the figure, it was found that a = —1.006(2) 
and z = —0.990(4). Given a and P the crossover expo- 
nent can also be obtained by solving z = a//3, yielding 
z = —1.001(4). Such result indeed agrees with our nu- 
merical result. Finally, with the values of the critical 
exponents, the scaling hypothesis can be verified. Fig. 
[3] (b) shows the collapse onto a single universal plot of 
five different curves of the P for different values of the 
dissipation parameter. 

We have also considered the behavior of the average 
standard deviation of /, which is defined as 



1 ^ /= 



(14) 
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(11) FIG. 7: (Color online) Plot of 2/ as a function of x for /3 = 0.5. 



We have considered an ensemble of M = 2000 differ- 
ent initial conditions iterated up to 10^. The variable /q 
has been fixed as /q = 10 -^iiT and ^0 were uniformly dis- 
tributed along ^0 ^ [0, 27r]. Figure[5](a) shows the behav- 
ior of the as a function of the number of collisions n for 
different values of the dissipation parameter and in Fig. [5] 
(b) their collapse onto an universal plot. The acceleration 
exponent is /3 = 0.501(3). In such a case ujgat oc 7^ and 
Tlx oc 7^, and cr = -0.489(1) and ^ = -0.970(3) as can 
be seen in Fig. [6l A similar scaling analysis shows a re- 
lationship ^ = cr//3 = —0.976(7), well in agreement with 
the numerical result for ^. This set of critical exponents 
is the same as observed in a dissipative bouncer model 
[29] where a classical particle collides inelastically with 
a time periodically moving wall under the effect of grav- 
ity. It also is the same set as observed for the oval-like 
billiard moving under the effect of a quadratic frictional 
force [30]. Although both, the systems and the kind of 
dissipation, are different, the phase transition they are 
experiencing, namely, suppressing the unlimited growth 
of the dynamical observable (energy for the oval-like bil- 
liard, deviation of the average velocity for the dissipative 
bouncer and deviation of the average action on the stan- 
dard), belongs to the same class of universality. 

The behavior observed in Fig. [3] (b) and Fig [5] (b) to- 
gether with the scaling discussion leads to an universal 
plot. This kind of behavior was also observed in many 
different systems, under different approaches for the con- 
servative and dissipative dynamics and considering differ- 
ent degrees of freedom mappings. To illustrate few ap- 
plications more than those discussed above in dissipative 
systems we nominate: (i) oval billiard [3l|; (ii) Lorentz 
gas [32] and; (iii) elliptical billiard [33^. For these cases, 
the transformation uo ujjujsat — V (or V V jVsat = V 
or P P/Psat) and n ^ n/ux = x dl\ the curves 
for all these systems collapse onto a single universal plot 
as can be seen in Fig. [3 (for j3 — 1/2) and Fig. [8] (for 
/3 = 1). Based on this, we propose the following empirical 
function 
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FIG. 8: (Color online) Plot of ^ = 1^ jl^sat as a function of 
X — nlux for the standard map where ^ —\. 



where the only parameter that remains is the acceleration 
exponent /3. As one can see in Fig. [3 and Fig. [8] the 
agreement of our function j{x) with the numerical data 
is quite good over many orders of magnitude with only 
some disagreement near the crossover. As mentioned, we 
have considered two values for the acceleration exponent, 
/3 = 0.5 and /3 = 1, which includes a number of different 
physical systems. 

III. CONCLUSION 

As the concluding remark, some results for a dissipa- 
tive standard map have been addressed. For the regime 



of strong dissipation, the model exhibits a period dou- 
bling bifurcation cascade, where the so called Feigen- 
baum's ^ = 4.6691(1) was numerically obtained. When 
weak dissipation is taken into account, for small values 
of the mixed structure of the phase space is changed, 
the elliptic fixed points are replaced by attracting fixed 
points. No large chaotic attr actors are observed in such a 
regime. For large the behavior of the average energy 
is considered in the framework of scaling. Once the expo- 
nents /3, OL and z are obtained, the scaling hypotheses are 
confirmed by the perfect collapse of all curves onto a sin- 
gle universal plot. The scaling behavior is also described 
in terms of the average standard deviation of the action 
variable. Given the exponents /3, a and ^ we can confirm 
that the dissipative standard map belongs to the same 
class of universality of a dissipative time dependent ellip- 
tic billiard (2-D) [33|, non-dissipative Fermi-Ulam model 
(1-D) [34], the corrugated wave guide (1-D) |35i], for the 
range of the control parameters considered. Finally, we 
propose an empirical function f{x) = x^/(l + x)^ in or- 
der to describe the universal scaling behavior observed 
when dissipation is taken into account for many different 
systems. This universal function gives a good descrip- 
tion of the empirical numerical data over many orders of 
magnitude. Using it, only one empirical parameter is left 
in the system, namely the acceleration exponent /3. 
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